Vortex structures of rotating spin-orbit coupled Bose-Einstein condensates 
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We consider the quasi-2D two- component Bose-Einstein condensates with Rashba spin-orbit (SO) 
coupling in a rotating trap. The rotation angular velocity couples to the mechanical angular mo- 
mentum which contains a non-canonical part arising from SO coupling. The effects of an external 
Zeeman term favoring spin polarization along the radial direction is also considered, which has the 
same form as the non-canonical part of the mechanical angular momentum. The rotating con- 
densate exhibits a variety of rich structures as varying the strengths of the trapping potential and 
interaction. With a strong trapping potential, the condensate exhibits a half-quantum vortex-lattice 
configuration. Such a configuration is driven to the normal one by introducing the external radial 
Zeeman field. In the case of a weak trap potential, the condensate exhibits a multi-domain pattern 
of plane- wave states under the external radial Zeeman field. 
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I. INTRODUCTION 

Spin-orbit (SO) coupling plays an important role in 
various aspects in condensed matter systems including 
spintronics [1[ and topological insulators 0, Q. How- 
ever, SO effects in bosonic systems has not been attracted 
much attention until recently. For example, 4 He atoms 
are spinless and ultracold spinful bosons are too heavy 
to exhibit relativistic SO coupling. This situation is sig- 
nificantly changed by the recent experimental progress in 
both semiconductor exciton systems and cold atom sys- 
tems with synthetic gauge fields. Excitons are composite 
bosons of electrons and holes. Their effective masses are 
light enough to exhibit relativistic SO coupling. Exotic 
SO coupled condensates with stripe and skyrmion types 
spin texture configurations was theoretically predicted by 
Wu and Mondragon-Shem [4]. Excitingly, spin textures 
have been observed in the SO coupled exciton conden- 
sates by High et al @. On the other hand, many theo- 
retical schemes have been proposed in ultracold atomic 
systems to create artificial non-Abelian gauge fields by 
using; laser-atom interactions which generate ef- 

fective SO coupling without special relativity. 

It has been shown that bosons with SO coupling sup- 
port exotic ground states beyond the "no- node" theo- 
rem jl8l-[2lj|. This theorem states that the ground state 
wavefunctions of bosons under very general conditions 
are positive definite, which is essentially a direct re- 
sult of the Perron-Frobenius theorem of matrix analy- 
sis (22[. However, the linear coupling to momentum in 
the SO coupling invalidates the proof of the "no-node" 
theorem. For example, spontaneous time-reveal symme- 
try breaking states exhibiting spin-density wave order- 
ing 0, [l9|, [28l43lj and spontaneous half-quantum vortex 
configuration 0, Gil have been studied. Both of them 
exhibit either nodal or complex- valued condensate wave- 
functions, and thus are beyond the "no- node" theorem. 
Especially, the realization of SO coupled Bose-Einstein 
condensations (BEC) of 87 Rb [32|,[33[ provides a valuable 



opportunity to investigate this type of exotic physics, ex- 
perimentally. Another way to bypass "no-node" theorem 
is to employ the meta- stable excited states, in which "no- 
node" theorem does not apply either. For example, cold 
alkali bosons have been pumped into the high orbitals in 
optical lattices [34|, [35| . It was shown that interactions 
among p-orbital bosons obey an "orbital Hund's rule" , 
which generates a class of orbital superfluid states with 
complex-valued wave functions breaking TR symmetry 
spontaneously [23l-[27|]. 

On the other hand, vortex properties in rotating BECs 
are a characteristic topological feature of superfluidity 
including 4 He and ultra-cold bosons, which have been 
studied extensively both experimentally and theoretically 
[3^ ]. For spinor BECs and spinful Cooper pairing super- 
fluidity (e.g. superfluid 3 He A and B-phases), exotic spin 
textures and fractional quantized vortices can form un- 
der rotation (38[. However, to our knowledge, the vortex 
properties of rotation SO coupled BECs have not been 
thoroughly investigated before. 

In this article, we investigate the rotating SO coupled 
condensate in a quasi-2D harmonic trap with the angular 
velocity along the z-axis. The angular velocity couples to 
the mechanical angular momentum whose non-canonical 
part behaves like a Zeeman term polarizing spin in the ra- 
dial direction. We also consider the effect from an exter- 
nal Zeeman term with the same form. The single particle 
ground states in the absence of interaction can have non- 
zero vortex numbers, which differ by one in the spin- up 
and down components as a result of SO coupling. With 
many-body interactions, the rotating condensate exhibit 
a variety of configurations depending on the strengths of 
the trapping potential and interaction. If the trapping 
potential is strong and interaction is relatively weak, a 
half-quantum vortex lattice is formed under rotation. Its 
spin configuration is a lattice of skyrmions. The conden- 
sate of the spin up component breaks into disconnected 
density peaks, which overlap the vortex cores of the spin- 
down condensate. The presence of the external Zeeman 
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field drives the system from a half quantum vortex lat- 
tice state to a normal quantum vortex lattice state. In 
the case of a weak trap potential, the condensate favors 
a plane-wave state or a two-plane-wave state with twist 
phase profiles under rotation. With the external Zeeman 
field, the condensate develops multi-domain configura- 
tion of plane- wave states. The configuration of wavevec- 
tors can be clockwise or counter-clockwise depending on 
the direction of the field. 

The rest part of the paper is organized as follows. The 
model Hamiltonian of the rotating Rashba coupled BEC 
is introduced in Sect. [TTJ The solution of the single parti- 
cle wavefunction is presented in Sect. IIIII The rich struc- 
tures of the vortex configurations with spin textures are 
given in Sect. [Vl onclusions are given in Sect. |Vl 



II. THE MODEL HAMILTONIAN 

We consider the quasi-2D two-component BECs with 
Rashba SO coupling in the xy-pl&ne subject to a rotation 
angular velocity Vt z along the z-direction. The free part 
of the Hamiltonian of Rashba SO coupling under rotation 
is defined through the standard minimal coupling as 



non-canonical part in L mech plays a crucial role during 
the understanding the single-particle ground state prop- 
erties. 

For the later convenience, we also introduce an external 
spatially dependent Zeeman term as 



H B = 



ex,x cr x + B eXi y(iy)^ v {r), (5) 



where B ex (f) = (Box, Boy, 0) varies linearly in the xy- 
plane. Such a term can tune the strength of the non- 
canonical part of the mechanical momentum, which ren- 
ders the model adjustable in a wider range of the param- 
eter space. 

Many efforts have been made to implement the above 
Hamiltonian in ultra-cold atomic gases. Several schemes 
have been proposed to generate Rashba SO coupling 
0, [H, G3 with tunable SO combing strength. In par- 
ticular, proposals in Ref. [HI, [I?} have the advantage 
to overcome the drawback of the spontaneous emission 
in the tripod scheme. The spatially dependent Zeeman 
term Hb can be generated through coupling two spin 
components using two standing waves in the x and y di- 
rections with a phase difference of tt/2. The resulting 
Rabi coupling is written as 



Ho = [ d 3 r^l(r)[^-(-ihW + MXz x a - A) 2 - fi ~ & [s'm(k L x) + i sm(k L y)] ^\(r)^(r) 

J ZlVl 



h.c. 



(6) 



V ext (7)--Mn 2 z (x 2 + y 2 )]^Ar), 



(1) 



where a = a x x + a y y + a z z with a x ,y,z the usual Pauli 
matrices; A is the Rashba SO coupling strength with the 
unit of velocity; /i, v take values of t, | as pseudospin 
indices; A — (— MQ z y, MQ z x, 0) is the vector potential 
from Coriolis force; V ext (r) = \Mujt(x 2 + y 2 ) is the ex- 
ternal harmonic trapping potential; the last term in Eq. 
[1] is the centrifugal force. The interaction part H int is 
defined as 



Hjnt — 



d 3 ^t(fO^(r)^(r)^(0- ( 2 ) 



We assume the equal intra-component interactions as 
#tt = 9\X = 9i an d inter-component interaction = gc 
with c a constant coefficient. 

Due to the presence of SO coupling, Q z couples to the 
mechanical angular momentum L mech rather than the 
canonical one L z . We extract this coupling from Eq. [T] 
as 



In the region of x, y <C 27r//cL, it reduces to the desired 
form of Eq. [5] with Bo — flkL- 



III. THE SINGLE PARTICLE SPECTRA 

We start with the non-interacting Hamiltonian Ho + 
Hb, to gain some intuition. The con fining trap is char- 
acterized by the length scale I = ^/h/Muj. We define 
another length scale l so = H/(MX) from SO coupling. 
The ratio between them a = l/l so is dimensionless pa- 
rameter to describe the strength of SO coupling. For the 
typical setup used in the NIST group [33|, a ~ 10. Below 
we vary the values of a from ~ 10. Experimentally, the 
regime of small a can be reached by using a deeper trap 
potential. 

If without the confining potential and rotation, the 
single-particle eigenstates is of the form 



Ak-r 



±,k 



±,k), 



(7) 



H^ot — 



J d 3 r^(r) 



\1V 



where 



jjmech 



MX(xa x +y<r y ). 



(3) 



(4) 



Therefore, rotation in the presence of SO coupling in- 
duces an effective magnetic field distribution Br(v) = 
Vt z MX(x, y, 0) in the x?/-plane. As we will see below, this 



where |=L,fc) = (l,^e z6> fc) T , and 0^ is the azimuthal 

angle of k. Since the condensate is uniform along the z 
direction, we always have k z = for the ground state. 
The corresponding dispersion relations come into two 
branches e± = h 2 (k 2 ± 2k k)/(2M) with k = l/l so - 
Therefore, the single particle ground states are infinitely 
degenerate along a ring in momentum space with radius 
k . 

The external harmonic potential has an important ef- 
fect which lifts the degeneracy along the Rashba ring as 
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FIG. 1: The canonical angular momenta m of the single par- 
ticle ground states described in Eq. [TUl'u.s. p for 7 = —0.1, 
0.0, and 0.15, respectively. 

pointed out in Ref. [4] . In the momentum representation, 
the harmonic potential becomes ^Moo 2 (ihV g — A') 2 in 
the lower branch and couples different plane wave states 
around the Rashba ring, where A'(k) = j*\Vz\ip_ r) 
corresponding to a 7r-flux at the origin. Therefore the 
motion along the Rashba ring is quantized and main- 
tains time-reversal (TR) invariance. The single particle 
spectra exhibit the fermion-type Kramer degeneracy with 
T 2 = —1. The lowest single particle eigenstates carry 
j z = =b|. As shown in Ref. Q, the angular quantization 
gives rise to the dispersion on j z as 

\\jz\ 2 tu*r. (8) 
or 

On the other hand, the radial quantization is the same 
as in the ordinary harmonic trap, which is at the order of 
huJT [4]. In the strong SO coupling limit, i.e., a ^> 1, the 
dispersion over j z is nearly flat. Thus the radial quantum 
number can be viewed as band index, and the quantum 
number j z marks each state in the band. 

To be more precise, we define two independent annihi- 
lation operators as a d — \(z + 2d z ) and a g = \(z + 2d z ) 
where z = (x±iy)/l and z is the complex conjugate of z 
[37l . [39j . The single-particle Hamiltonian can be rewrit- 
ten in the unit hu as 

H + H B = (1 - p)N d + (1 + p)N g + 1 

+ a {[(1 - K)a d - (1 + ^)aj]cr + + h.c.} , (9) 

where 

p = il z /u, N d = ci ] d a d , N g = a\a g , a + = ^(a x + zcr y ), 

and k = j + p with 7 = Bq/(Muj\). The corresponding 
canonical angular momentum reads L 2 = ftZ^ = ^(7V^ — 
N g ). The A^-term represents the combined effect from the 
non-canonical part of H rot and the Zeeman term He- 

We diagonalize Eq. [9]to obtain the single particle spec- 
tra, and present the solutions in the coordinate represen- 




FIG. 2: Density profiles of spin- up and down components 
for the single-particle ground state 6- _ 1 with the parameter 

Jz— 2 

values of a = 4, and p — 7 = 0. Only the spin-down compo- 
nent carries a vortex. The corresponding density profile for 
4>j z= _i is obtained by interchanging the spin indices. 

tation, in which the ground state wavefunction reads as 

The total canonical angular momentum j z = l z + \o z = 
m + I remains a conserved quantity, thus the canoni- 
cal orbital angular momenta in the two spin components 
differ by one due to SO coupling. Fig. [1] shows m as a 
function of the rotational angular velocity p for different 
external magnetic field B ex at a = 4. In the absence 
of B exi the total angular momentum j z = — ^ for small 
p and decreases when p — » 1. Introducing the field B ex 
changes the ground state dramatically. If B ex is par- 
allel to the induced magnetic field Br, i.e., 7 > 0, j z 
first decreases then increases with the rotational angu- 
lar velocity p. However, for 7 < 0, j z increases with p 
monotonically. 

The above results can be understood as follows. In the 
case of Q z = 0, the two states </>^ =± i are degenerated 
due to TR symmetry. Since only one of the two spin com- 
ponents carries a vortex, the ground state can be viewed 
as a half quantum vortex state with the density profiles 
of two spin components shown in Fig. O The spin den- 
sity distributions exhibit skyrmion-type texture configu- 
rations, as depicted in Fig. [3] A and B. Intuitively, one 
might expect that an infinitesimal Q z selects the 4 > j z = ± 
state since it has lower rotational energy —Q Z (L Z ). How- 
ever, the presence of the induced magnetic field Br con- 
tributes another term to the total rotational energy of 
the system as 

(H rot ) = -n z (L z )-B R .(a). (11) 

The spin pattern (a) for <j>j z= * in the x?/-plane is anti- 
parallel to Br near the trap center (see Fig. EfA)), which 
is energetically unfavorable. Therefore, when —Br • (a) 
dominates, j z of the ground state can be — \ for a rotating 
trap. As increasing Q z , the condensates expand, which 
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FIG. 3: The spin density vector distributions along the x- 
axis lie in the xz-plane as shown in (A) j z = | and (B) 
j z — — \ with a = 4, and p = 7 = 0. They are time-reversal 
counterparts to each other, and both exhibit the skyrmion- 
type texture configuration. 



also favors the magnetic energy term. The total angular 
momentum j z can decrease even when p increases. Such 
a counter-intuitive effect for the ground state constitutes 
a characteristic feature of SO coupled BECs in a rotat- 
ing trap. Introducing the external magnetic field B ex 
strengthen or weaken this effect induced by Br depend- 
ing on its direction, which explains the different behaviors 
of m with p for 7 > and 7 < 0, as shown in Fig. [TJ 



IV. VORTEX CONFIGURATIONS OF 
ROTATING SO COUPLED BEC 

Interaction effects in the absence of rotation have been 
investigated extensively in the literature, which are sum- 
marized below. In the case of a strong trapping po- 
tential and weak interaction, the single-particle energy 
dominates. The condensate maintains rotational symme- 
try but spontaneous breaks TR symmetry 0, [l9j . One 
spin-component carry one vortex, and the other is non- 
rotating, thus the condensate possesses a half-quantum 
vortex. The total angular momentum of each particle 
is \j z \ — \- I n momentum space, this kind of ground 
state distributes uniformly around the Rashba ring. On 
the contrary, if the trapping potential is weak and inter- 
action is strong, the condensate breaks rotational sym- 
metry. The condensate is approximately superposition 
of plane-wave states modified by the cylindrical bound- 
ary condition. Results based on the Gross-Pitaevskii (G- 
P) equation show that the spin-spiral condensate with 
two counter-propagating plane- waves is favored at c > 1 , 
while a single plane- wave is favored at c < 1 (281431). 
These two different condensates are degenerate for the 
spin-independent interactions, i.e., c = 1. However, 



calculations including quantum fluctuations of the zero- 
point energy show that the spin-spiral state wins at c = 1, 
and thus shift the phase boundary to a smaller value of 

In this section, we study the vortex configurations of 
SO coupled BECs in both cases. The results of strong 
trapping potentials and weak interactions are presented 
in Sect. IIV A[ and those of the opposite limit are pre- 
sented in Sect. HVBl 

A. Vortex lattice configurations with a strong 
trapping potential 

In this subsection, we turn on rotation and consider a 
strong trapping potential with a small value of a. The 
ground state condensate is obtained by numerically solv- 
ing the following SO coupled GP equation. We assume 
that the condensate is uniform along the z-axis, and de- 
fine the normalized condensate wavefunction (^,^) T 
satisfying J d 2 f{\^\ 2 + \^\\ 2 ) = 1. The dimensionless 
version of the G-P equation can then be written as 

= f^ + f3(\^\ 2 +c\^\ 2 )^, (12a) 
£^ = T^ v + /?(|^| 2 + <# t | 2 )^, (12b) 

where 

f = -h 2 (dl + dl) + al(-id y a x + id x G y ) 

+ ^(x 2 + y 2 ) - p(~ixdy + iyd x ) 

- —(xcr x + y(j y ), (13) 

where p is the chemical potential; the interaction param- 
eter f3 = gN/(hul z ); N is the particle number in the 
condensate; l z is the system size along the z-direction. 

The density and phase configurations at various pa- 
rameters are shown in Fig. |4](a) - (g), which exhibit rich 
structures of vortex-lattice. We look at Fig. 3] (d) in 
the absence of B eXl i.e. 7 = 0. The density distribution 
of the spin-up component is composed of several discon- 
nected density peaks near the trap center. On the other 
hand, the low density region is connected in contrast to 
the usual vortex lattice structure in which the low density 
region of vortex cores is disconnected. Nevertheless, we 
identify the locations of the singular points of the phase 
distribution pattern around which the phase winds with 
an integer number. These singular points are squeezed 
out to the edge of the condensate. On the other hand, the 
spin-down component exhibits the regular vortex-lattice 
structure, whose vortex cores overlap with the density 
peaks of the spin- up component. Around each vortex 
core, the two spin components show a half-quantum vor- 
tex configuration as those depicted in Fig. [2j Therefore, 
the condensates of two components together exhibit a 
lattice of half-quantum vortices. The corresponding spin 
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FIG. 4: From left to right: the density and phase profiles 
of spin-up and down components with parameter values of 
a — 0.5, /3 = 10, p — 0.97, and c = 1. From (a)-(g), 7 is taken 
as 0.5, 0.25, 0.1, 0.0, -0.1, -0.25, and -0.5, respectively. 
At small values of I7I in (c) - (e), a half quantum vortex 
lattice is formed near the trap center. The spin-up component 
breaks into several density peaks, and the low density region 
is connected. As increasing the magnitude of |^y| (b) and (f), 
the half-quantum vortex lattice evolves to the normal vortex 
lattice. For the large value of | = 0.5 (a) and (g), the 
condensates show a lattice configuration around a ring. The 
black circle with an arrow indicates the the direction of the 
circulation around the vortex core. 



density vector (a) shows a skyrmion-lattice structure, as 
shown in Fig. [5] 

Now we turn on the external Zeeman term Eq. [5] For 
both cases of 7 > and 7 < 0, at small values of |7|, 
the half-quantum vortex lattice still forms, which is sim- 
ilar to that at 7 = as depicted in Fig. 2] (b, c, e). 
As increasing the strength of B exi i.e., |7|, more vortices 
appear as depicted in Fig. 2] (b,f). The condensates 
of the spin-up component gradually evolves to the usual 
vortex-lattice configuration. The high density region be- 
comes connected, while the density minima become dis- 
connected vortex cores. On the other hand, the con- 




FIG. 5: The ground- state spin density vector of Fig. [4] (d) 
with parameter values of a — 0.5, f3 = 10, p = 0.97, c = 1, 
and 7 = 0. The projection of (a) in the xy-plane is shown 
as black vectors. A color map is used to illustrate the (a z ) 
component. 



densates of the spin-down component remains the usual 
vortex lattice configuration. For even larger values of |7|, 
the ring-shaped vortex lattice with a giant vortex core is 
observed as shown in Fig. [4] (a) and (g). This is be- 
cause the combined effect of the harmonic trap V ex t(r) 
and the additional Zeeman term Hb shifts the poten- 
tial minimum to a ring in real space with the radius of 
r = ajl = \Bo\/(Muj 2 ). The condensates of both spin 
up and down components distribute around this ring and 
from a giant vortex configuration. Additionally, the Zee- 
man term grows linearly as increasing r and favors in- 
plane polarization of S. As a result, the vortex cores 
of the spin up and down components overlap with each 
other. 

We stress that in all cases in Fig. [H (a-g), the vor- 
tex numbers in the spin-up and down components differ 
by one, which is a characteristic feature brought by SO 
coupling. As shown in Eq. [lOl for the eigenstate of the 
single-particle Hamiltonian with j z = m + |, the two 
spin components carry different canonical orbital angu- 
lar momenta m and ra+1, respectively. In the presence of 
interaction, the giant vortex splits into a lattice of single- 
quantum vortices in each spin component. Nevertheless, 
the total vortex number in each component remains un- 
changed and differs by one. 



B. Weak trapping potential 

In this subsection, we study the rotating SO coupled 
BEC with a weak trapping potential and strong interac- 
tions. 

Fig. [6] shows the density and phase profiles of each 
spin component in the absence of external magnetic field 
B, i.e., 7 = 0. In Fig. [6] (a) with c < 1, the conden- 
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FIG. 6: From left to right: the density and phase profiles of 
spin up and down components with the parameter values of 
a = 4, j3 = 20, and 7 = 0. (a) c = 0.6 and p — 0.1, a plane- 
wave-like state is obtained with a distorted phase pattern; (b) 
c = 1.2 and p = 0.1, the spin- spiral condensate is favored; (c) 
c — 1.2 and p = 0.5, the condensate exhibits an intermediate 
configuration between those of (a) and (b). The color scales 
for the density and phase distributions are the same as those 
in Fig. H 



sate is a twisted plane- wave state subject to the cylin- 
drical boundary condition. The spin polarization mainly 
lies in the xy-plane. In the representation eigen-basis of 
s z , the spin up and down components show nearly the 
same distributions of density and phase profiles. Never- 
theless, the phase distribution is distorted from the ex- 
act plane-wave state. On the other hand, as depicted in 
Fig. E](b), at c > 1 the spin-spiral-like condensate with 
two counter-propagating plane- waves is still favored with 
twisted phase profiles. As shown in Fig. [6](c), increas- 
ing the angular velocity p gives rise to an intermediate 
configuration between the distorted spin-spiral and the 
single-plane wave states. In all the patterns, vortices lo- 
cate either on the edge of the condensate or the density 
minima of each component. 

Next, we consider the case of 7 7^ 0. Introducing Hb 
significantly enriches the structures of the rotating SO 
coupled condensates. We only consider a small angu- 
lar velocity at p = 0.1 for the reason of numerical con- 
vergence, but vary the values of 7 from 0.5 ~ —0.7 as 
presented in Fig. [7] (a)-(h), respectively. With small 
and intermediate values of I7I (e.g. Fig. (b)-(f)), the 
condensate breaks into several domains. Inside each do- 
main, the condensate can be approximated as a single 
plane- wave state. Vortices center around the local den- 
sity minima. The local wavevectors are configured such 
that the local spin polarization (S) align along the local 
Zeeman field of B ex (f). If 7 > at which the external 
Zeeman field enhances the rotation induced ones, we ob- 
tain a clockwise configuration of wavevectors. There is 
one more vortex with the negative phase winding in the 
spin up component than in the spin down component, 
which reflects the "anti-paramagnetic" feature. On the 
contrary, if 7 < 0, the anti-clockwise patterns of wavevec- 
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FIG. 7: From left to right: the density and phase profiles 
of spin up and down components with parameter values of 
a = 4, f3 = 20, c = 1, and p — 0.1. From (a)-(h), 7 is taken 
as 0.5, 0.3, 0.1, -0.05, -0.25, -0.35, -0.6, and -0.7, respec- 
tively. The black arrow in each domain represents the local 
wavevector direction of the corresponding plane- wave state, 
which shows a clockwise or counter-clockwise configuration 
depending on the sign of 7. For sufficiently large values of 
I7I, condensates distribute around a ring in space forming a 
giant vortex. The color scales for the density and phase dis- 
tributions are the same as that in Fig. [4] The black circle 
with an arrow indicates the the direction of the circulation 
around the vortex core. 
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tors is favored. Similarly, the spin-down component also 
carries one more vortex than the up component. 

At small values of I7I, two domains are formed as de- 
picted in Fig. [3(c) and (d). The vortices organize into 
straight-lines between two domains. A variational wave- 
function is constructed as 



e -r 2 /(2a 2 ) 



(14) 



where without loss of generality, we choose the wavevec- 
tor /co = koe y ; a is radius of the condensate; 6 is the 
relative phase difference between the two plane wave do- 
mains; |/_ ?+ (x)| 2 = (e± x l w -f are smeared step 
functions with W the width of the domain wall. We 
assume a ^> (W, 1/fco)- Such variational wavefunction 
has neglectable contribution to the energy term (H rot ). 
This explains why the two domain pattern is absent by 
increasing the rotational angular velocity p only, but ap- 
pears immediately even at small values of |7|. With in- 
creasing I7I, the condensate breaks into more and more 
domains as in Fig. 0(b), (e) and (f). 

As further increasing I7I, domains connect together as 
a giant vortex as shown in Fig. [7] (a, g, h). The conden- 
sates of both spin up and down components distribute 
around a ring with the radius of a\j\l and overlap each 
other. This is a giant vortex configuration with a texture 
of spin aligned along the radial direction. The phase 
winding numbers of the spin-up and down components 
differ by one due to the SO coupling. 



V. CONCLUSION 

To summarize, we have considered the vortex struc- 
tures of SO coupled BECs in a rotating trap combined 



with an external spatially dependent Zeeman field. In 
the case of strong confining potentials and weak inter- 
actions, the condensate exhibit vortex-lattice structures. 
As varying the magnitude of the external Zeeman field, 
the configuration evolves from a half-quantum vortex- 
lattice to a normal one. In the opposite limit, the conden- 
sate develops multi-domain patterns with the external 
Zeeman field. Each domain represents a local plane- wave 
state, whose wavevector exhibit a clockwise or counter- 
clockwise configuration. Domain boundaries play the role 
of like vortices. 
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Note added Near the completion of this manuscript, 
we notice a recent paper studying the rotating Rashba 
SO coupled BEC, which considered a special case in the 
presence of the extra term of Eq. [5] with 7 = — p [40] . Our 
work has studied the general cases, including the pure 
rotation without the external fields which corresponds to 
7 = 0. 
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